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We introduce a general method for designing taifored lattices of magnetic microtraps for ultracold 
atoms, on the basis of patterned permanently magnetized films. A fast numerical algorithm is used 
to automatically generate patterns which provide optimal atom confinement while respecting desired 
lattice symmetries and trap parameters. The algorithm can produce finite and infinite lattices of 
any plane symmetry; we focus specifically on square and triangular lattices which are of interest 
to future experiments. Typical trap parameters are discussed as well as the impact of realistic 
imperfections such as finite lithographic resolution and magnetic inhomogeneity. The presented 
designer lattices open new avenues for quantum simulation and quantum information processing 
with ultracold atoms on atom chips. 
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I. INTRODUCTION 

Ultracold atoms in periodic potentials are ideal sys- 
tems for exploring a diverse range of new physics, from 
the behavior of individual trapped atoms to rich many- 
body phenomena. At the core of most current exper- 
iments are optical lattices, produced by the interfer- 
ence of intersecting laser beams [IHS], which have pro- 
vided for example unprecedented access to Josephson ef- 
fects squeezed states [S], low-dimensional quantum 
gases [1] and novel quantum phases [J- Particularly ex- 
citing applications arise in the context of quantum sim- 
ulation fE^ and quantum information processing |S1 llOj , 
where atoms or molecules trapped in periodic arrays hold 
great promise for the storage and manipulation of quan- 
tum information. Lattices of magnetically trapped atoms 
or ensembles of atoms, interacting via a switchable Ry- 
dberg excitation, are promising systems for the gener- 
ation of cluster states and other many-particle entan- 
gled states. Thus they could provide a suitable platform 
for implementing measurement-based quantum informa- 
tion processing schemes as well as exploring fault tolerant 
schemes [TT^14) . 

For practical devices, an attractive alternative to op- 
tical lattices are arrays of magnetic microtraps produced 
by current-carrying wires or patterned magnetic films 
on atom chips |15H21j . In principle, these magnetic lat- 
tices provide robust potentials for manipulating atoms, 
combined with a high degree of simplicity and excellent 
design flexibility. For example, inter-atomic distances 
are not restricted to certain fractions of optical wave- 
lengths, allowing for a trade-off between inter-site cou- 
pling strength and single-site resolution. However, de- 



spite these advantages, only one-dimensional and a few 
two-dimensional lattices have been successfully imple- 
mented or proposed to date. The problem of finding the 
necessary magnetization pattern, or wire shape, which 
generates a desired set of atom traps currently relies on 
experience and trial-and error; what has been lacking is 
a general procedure for designing lattices with arbitrary 
arrangements of traps with the desired properties 

In this paper we present a general method for creating 
tailored lattices of magnetic microtraps by controlling the 
geometric patterns of perpendicularly magnetized planar 
films. We employ a linear programming algorithm to 
find optimal single-layer magnetization patterns which 
produce desired lattice symmetries with specified trap 
parameters. This algorithm is similar to one previously 
used for the optimization of surface-electrode ion-trap 
lattices based on radio- frequency (rf) electric fields |22) . 
It can be used to create a wide variety of designer-lattice 
geometries with arbitrary trap arrangements, opening 
new avenues for simulating condensed matter systems 
and for quantum information processing with ultracold 
atoms on atom chips. The planar patterns generated by 
the algorithm are simple to fabricate and can be straight- 
forwardly implemented in existing experiments. Current- 
carrying wires tracing the edges of the optimized patterns 
represent an alternative experimental implementation of 
our results [23]. 

In the first part of the paper we give an overview of 
the defining properties of magnetic films (Section |II A ) 
and lofFe-Prit chard traps (Section II B), followed by a 



description of the linear programming algorithm (Sec- 
II C|). We then turn to the results of the algorithm 



tion 



nil. To highlight the flexibility of the method 
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(Section 

we focus on two desirable lattices we can generate 
square lattice and one triangular lattice, of interest for fu- 
ture experiments. These lattices offer tight confinement 
and a high degree of symmetry. In Section |IV| we con- 
sider typical trap parameters produced by these lattices 
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and account for realistic imperfections such as finite litho- 
graphic resolution and inhomogeneity which could arise 
during fabrication. Section [TV D| discusses how these mi- 
crotrap lattices can be loaded. 



II. THEORETICAL FRAME 

A. Permanent magnetic microstructures 

For our present purposes of designing magnetic po- 
tentials we consider single-layer patterned magnetic thin 
films with perpendicular magnetization |19| 121) . There is 
an analogy between such patterned magnetic films and 
an effective boundary current which would produce an 
equivalent magnetic field [53^. We describe magnetic 
fields as gradients B{r) — — V^(f) of the scalar magnetic 
potential \E'(r) [23]. In particular, the scalar magnetic po- 
tential generated above a thin (typical thickness S <^ z) 
permanently out-of-plane magnetized layer of material in 
the xy plane V is 

l-^ioSM, [ m{x',y')Gix-x',y-y',z)dx'dy', (1) 



where /ip is the permeability of free space and is the 
remanent out-of-planc magnetization. Typical magnetic 
films achieve permanent out-of-plane magnetizations of 
the order of ~ 10^ A/m, and typical film thicknesses 
are of the order of (5 ^ 100 nm |19| (giving a magnetiza- 
tion current, or edge current in a wire implementation, of 
6Mz « 0.1 A); typical coercivities are around 500-1000 G, 
such that any external fields applied to create the mi- 
crotraps do not cause remagnetization. The dimension- 
less function m(x, y) represents the spatial dependence of 
the magnetization current (i.e., the product of remanent 
magnetization and film thickness); we assume it to vary 
between and 1 without loss of generality. The Green's 
function is 



G{x,y,z) = 



2tt{x^ + y2 + z2)3/2 



(2) 



and serves to propagate the planar Dirichlet boundary 
condition \E'(x,?/,0) = \iio5Mzm{x,y) through space 
while satisfying Laplace's equation V^^(r) = 0. 

We restrict our analysis to single-layer magnetization 
patterns with binary "step-like" magnetization (or thick- 
ness) variations, as they are relevant for current experi- 
ments due to the ease of patterning We note that 
the form of Eq. ([T]) is proportional to that which prop- 
agates an electric potential from an electrode plane into 
space [23]. Further, the shape constraints on a magne- 
tization pattern are mathematically equivalent to those 
on rf electrodes for ion trapping, as will be shown in Sec- 
tion [lIBj This leads us to the conclusion that the opti- 
mization algorithm of Ref. [21] , recently used to generate 
electrode patterns for ion trapping, is also ideally suited 



to designing new geometries of magnetic lattices for ul- 
tracold neutral atoms. In this sense the main goal of 
the remainder of this article is to find specific magnetiza- 
tion patterns m{x, y) which produce magnetic lattices of 
optimal strength (i.e., the stiffcst and deepest magnetic 
traps for fixed magnetization current SMz), given a set 
of constraints on trap arrangement and geometry. 

Despite the mathematical analogies, there are several 
practical differences between the designs of rf ion traps 
and magnetic microtraps. These differences concern the 
spatial derivatives of the electric or magnetic scalar po- 
tentials necessary to produce stable trapping potentials. 
In what follows we detail the conditions necessary for 
generating magnetic trapping geometries of interest to 
current experiments. 



B. lofFe— Pritchard traps 

Magnetostatic traps confine ultracold atoms of a spe- 
cific hyperfine state \F,mp) through a spatially varying 
Zeeman shift pseudo-potential |25j 



V{r) = mFgFf^B\\B{f)\\, 



(3) 



where (?f is the level's Lande factor and /ie is the Bohr 
magneton. Unlike for rf ion traps , here the trapping 
minimum of V{r) must be nonzero in order to prevent 
Majorana spin flips [IS]. loffe-Pritchard (IP) type traps 
are the simplest traps which guarantee this property. 

Let us assume that we wish to trap atoms at several 
points 7^^^ in space, with i = 1,2,3,... indexing the var- 
ious potential minima. Since the properties of each atom 
trap depend mostly on the local shape of the magnetic 
field, we require precise ways of specifying this shape in 
order to generate a desired lattice of atom traps. We base 
our analysis on a series expansion of the scalar magnetic 
potential around the trap points. 



w)-<''+j:-^\n~rr') 
i=i 

+ 1 t <l,.ir.-^)ir,-^)ir,-^) + 



(4) 



where (ri,r2,r^) = {x,y,z) are the Cartesian coordi- 
nates. The vectors -u^^-* have components uf\ while u^^-* 
and lo^^-* are tensors with components v^^j and w^fj f., re- 
spectively. In this notation, u^^^) = —B{r^^'>) specifies the 
local magnetic field at the desired trap position (apart 
from external bias fields), while v^^'^ is the local gradient 
tensor and w''^^ specifies the local magnetic field curva- 
tures. As described in Ref. J25j the coefficients of Eq. Q 
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must be fully symmetric under permutation of indices 
and satisfy 



^1,1 



-"2,2 



W 



i.1.1 



■ W, 



it) 
i.2,2 



W 



it) 
"3,3 

w 

4,3,3 



(5a) 
Vie {1,2,3} (5b) 



since V^^* = —V ■ B — 0. We use these vectors and 
tensors u^^\ v'^^\ and w^^^ as input parameters for our 



optimization algorithm (Section II C I in order to specify 



the shape and orientation of the magnetic field at the trap 
locations. In what follows we detail how these derivatives 
of * are constrained to specific forms for representing a 
particular desired lattice of IP traps. 

The defining characteristics of IP traps are that each 
gradient tensor v^-^^ has a zero eigenvalue associated with 
its 'loffe axis', and that the total local magnetic field 
('loffe field' Bi) is parallel to this axis: ■w(^) • Bi = 0. We 
therefore fix the field gradient tensors v'^^ at the trap 



locations [for concrete examples we use Eq. (13)]. 



Although it is possible to produce self-biased magnetic 
microtraps which require no external fields [26] . in prac- 
tice a homogeneous external bias field Bq is usually ap- 
plied to produce the loffe field Bi simultaneously at all 
trap sites. This dramatically simplifies loading the micro- 
traps (see Section IV D) and ensures that the magnetic 



potential remains finite away from the film surface (al- 
lowing for deeper confinement). To ensure all traps have 
equal depths for a given bias field Bq, we typically con- 
strain the first-derivative vectors m*-^' to be equal, and 
align the loffe axes at all trap positions r^^-'. However, 
there are no a priori constraints on the actual values of 
these vectors and tensors. 

If desired, additional constraints can be placed on the 
remaining two eigen-directions or the eigenvalues of the 
v^^\ on the local fields u^^\ or on the curvature ten- 
sors 10^^' in order to tailor the shapes, orientations and 
aspect ratios of the resulting trapping potentials. In par- 
ticular, the curvature matrix of an IP trap, which defines 
the characteristic oscillation frequencies (trap stiffness), 
is given by 



d^\\B{r)~B{f 



drjdri 



Bi 



I St 



(6) 

It is important that all eigenvalues of these curvature 
matrices be positive for a certain homogeneous bias field 
Bo = —B{f^^^) + Bi, in order to simultaneously trap 
atoms at all the positions r^^-' . Unlike for rf ion traps [52] , 
the three principal curvatures of IP traps are quite inde- 
pendent; prolate, spherical, oblate, and triaxial traps can 
be designed by judiciously choosing the conditions on the 
v^^^ and tt)^^) as well as operating at specific loffe field 
strengths. 

All constraints on the derivatives of are only speci- 
fied relative to each other, not by their absolute values. 
The optimization algorithm of Section |II C| aims at dis- 
covering the magnetization pattern which satisfies the 



constraints outlined in this section while maximizing the 
absolute values of the derivatives of VP. The fact that 
Eq. ^ depends on the loffe field strength makes it diffi- 
cult to introduce a dimensionless quantity describing the 
absolute strength of an optimized IP trap. However, it 
appears that jointly maximizing the v^^^ and to^^-' using 
our algorithm generally leads to deeper traps with higher 
trap frequencies, for fixed other physical constraints like 
the magnetization current SMz- 



C. Optimization algorithm 

The algorithm for finding the optimal magnetization 
pattern implementing a desired loffe-Pritchard trap ar- 
ray follows from Ref. [32]. We present it here in a more 
extended form. 

The first step in optimizing the magnetization pattern 
m(x, y) consists of subdividing the atom chip surface into 
N small domains. The shapes of these domains are ir- 
relevant for our purposes, and their number only serves 
to increase the resolution of the final magnetization pat- 
tern. It is not necessary that these domains cover the 
entire atom chip, that they lie in a single plane, or that 
they form a simply connected subset of the plane. For 
finite sets of traps the optimization may be performed on 
a finite-sized pattern in order to fully account for bound- 
ary effects. For large periodic lattices of traps it is useful 
to divide a single unit cell of the wallpaper group with 
the desired lattice symmetry into domains, and then as- 
sume that the same pattern is repeated indefinitely over 
the atom chip surface. 

We assume that within a domain a covering an area 
Va of the atom chip surface, the scaled magnetization 
m{x, y) — nia is constant. The total scalar magnetic 
potential, Eq. ([T]), can then be decomposed into a sum 
over domains. 



(7) 



where the functions 



i'aix,y,z) 



G{x-x',y-y\z)dx'dy' (8) 



are the scalar magnetic potentials induced by the do- 
mains indexed by a, assuming unit magnetization (since 
the term ^fioSMzfUa has been factored out). The present 
algorithm can also be applied to different problems such 
as non-planar geometries, electrode patterns, or hybrid 
atom-ion traps by replacing Eq. with more general 
expressions, or by evaluating Eq. Isl) numerically using 
finite-element or boundary-element methods ('S7] . The 
derivatives of the scalar magnetic potential are similarly 
parametrized 



<9«'(r) 



(9) 
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etc. With these definitions the hnear conditions of Sec- 



tion II B can aU be expressed in terms of the scaled do- 
main magnetizations vria without knowing their values 
a priori. We define m as the vector with components 
rria , which represents the optimal magnetization pattern 
to be found, restricted to TOq = or = 1 for non- 
magnetic or fully magnetized domains, respectively. Us- 
ing this notation, each linear condition of Section |IIB| 
can be formulated in the form Sk ■ rn = Cbk, where C is 
a common prefactor in all conditions and k indexes the 
different constraints. As in Ref. [55] the conditions of 
Section |IIB| are never imposed as absolute numbers but 
only relative to each other, while their absolute magni- 
tudes (through the prefactor C) are jointly maximized by 
the algorithm in order to achieve maximal strength of the 
Zeeman pseudo-potential for a given magnetization cur- 
rent SMz- All these conditions can be jointly expressed 
as a matrix equation 



A • m = Cb, 



(10) 



with \C\ to be maximized. 

In periodic arrays of microtraps the magnetization pat- 
tern is most naturally expressed in terms of its Fourier 
amplitudes rh{kx,ky) — T[m{x,y)]. Instead of using the 
Green's function of Eq. ([2]), the propagation of the scalar 
magnetic potential away from the atom chip then follows 
from the exponential damping of each Fourier compo- 
nent, 



^ m(fc„ fcj,)e'('=-"+'=''^)-"V^^ 



(11) 



where Re[. . .] refers to the real part. Since the Fourier 
transform is a linear operation, these Fourier amplitudes 
rh{kx, ky) can be expressed as linear functions of the do- 
main magnetizations TOq. Any linear constraints on the 
Fourier amplitudes can therefore also be brought into the 
form dk-rn — bk and thus included in Eq. ( 10 1. However, 



we will not be using such direct Fourier constraints in 
our examples (Section III I. 

The physical constraints describing the design for an 
atom trap lattice are thus formulated in the matrix equa- 
tion of Eq. (10 1. If the number of domains N is large. 



then this linear system of equations is highly underdeter- 
mined. Which of the solutions should we choose? The 
inhomogeneous solution mg = • b, calculated from 
the pseudo-inverse A^, incorporates the main structure 
of the magnetization pattern we will generate. However 
rriQ does not satisfy the constraint that rria € {0,1}: its 
components are not normalized in any useful way, and 
they take on a continuum of values. To remedy this we 
decompose the target solution into rh = Cmg + rh' with 
A-rn' = 0. Any vector of this form satisfies Eq. ( [To| ); but 
it serves to illustrate that the solution m should consist 



of the maximum possible contribution in the inhomoge- 
neous direction toq (since we wish to maximize |C|) plus 
just the right amounts of homogeneous solutions (in m') 
in order to satisfy the binary constraints. Now we de- 
fine A ~ A ■ (l — momJ/||TOo||^) = A — 6mJ/||mop, 
for which A-rn — A-rn' — because rn! - mo — 
(\A^Y ■ "^') • ^ = 0. The optimization problem can then 
be expressed as a homogeneous binary integer linear pro- 
gram (A, mo, {0, 1}), which schematically reads thus: 

(i) We seek a vector fn within the space satisfying the 
linear constraints A • m — 0. This gives a set of 
homogeneous conditions on m which do not involve 
the constant C, and which express the fact that the 
solution rh must satisfy the design constraints for 
the atom trap configuration, Eq. (10). 



(ii) All components of the vector m must satisfy ma S 
{0, 1}. This condition expresses the desire for a bi- 
nary magnetization pattern, consisting only of non- 
magnetic domains {ma = 0) and fully magnetized 
domains {ma = 1)- 

(iii) We wish to maximize \C\ = \fh - mo|/||mo|P, the 
"strength" of the solution in the inhomogeneous 
direction satisfying the design constraints. Larger 
values of |C| yield stronger atom confinement for a 
given magnetization current SM^. 

Such integer linear programs are known to be NP-hard 
to optimize |28j . if an exact solution even exists. For- 
tunately, and perhaps surprisingly, we can relax condi- 
tion ([ii| to < ma < 1 Va = 1 ... TV without loss of 
accuracy but drastically decreasing the computational 
complexity to 0{N), making large numbers of domains 
(high resolution) feasible. After relaxing condition ^ 
we are dealing with a linear program (A, mg, ... 1), for 
which very efficient algorithms are available [29] . The so- 
lution to the linear program is known to be globally opti- 
mal [2n]j and the obtained solutions consist of relatively 
large uniform patches of magnetization railed at either 
m{x, y) — or m(x, y) = 1. Further it is known that the 
number of domains a which are not railed is equal to the 
number of independent constraints in Eq. (10 1 ;22J. Since 



the number of domains can be made arbitrarily large, we 
are thus guaranteed to find solutions in which the ratio 
of un-railed to railed domains can be made vanishingly 
small. In practical cases the few un-railed domains usu- 
ally come to lie at the edges between large railed patches 
and can be safely rounded to or 1. 



1. limitations of the algorithm 

The optimization algorithm of the preceding section is 
a powerful new tool for designing magnetic lattice poten- 
tials with arbitrary geometries and trap characteristics. 
There are, however, some limitations which deserve men- 
tion. 
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Although any magnetization pattern generated by the 
above algorithm is guaranteed to produce the desired lat- 
tice of magnetic microtraps, it disregards the shape of the 
Zeeman pseudo-potential between microtraps. While this 
point is not critical in ion trapping applications [22 , it 
may be of great interest for atom trapping to have efh- 
cient tunneling paths between microtraps. To this end 
it may be necessary to add further linear constraints on 
the scalar magnetic potential. 

Generally, for trap stability the resultant Zeeman 
pseudo-potentials should not contain any points of zero 
magnetic field in the vicinity of the traps, since they will 
lead to trap losses due to tunneling from the IP traps and 
subsequent Majorana spin flips [25 . It is therefore more 
important than in ion trap setups that after a magneti- 
zation pattern m{x,y) is optimized, we check the corre- 
sponding Zeeman pseudo-potential 

Vir) = mFgFf,B\\B{r) - B(r<^)) + Bi\\ (12) 



2. implementation 

The linear optimization algorithm has been imple- 
mented independently in the Mathematica and Matlab 
programming environments for the general case of two- 
dimensional Bravais lattices, where each unit cell consists 
of ~ rii X n2 identical parallelogram-shaped domains. 
Further, a Mathematica implementation allows optimiz- 
ing non-periodic (finite) patterns with arbitrary polyg- 
onal domain shapes. All implementations use primal- 
dual interior-point algorithms for solving the linear pro- 
grams [53 and yield consistent results for the obtained 
optimal magnetization patterns. The typical number of 
domains (or Fourier modes) used in the periodic opti- 
mization is A^ ^ 200 X 200 which yields smooth mag- 
netization patterns and is a sufficiently high resolution 
to ensure that the corresponding magnetic potentials are 
precisely defined. 



for spurious null points. Such points can appear at unex- 
pected locations which depend on the external bias field 
Bo = -B{r^^'i) + Bi, and can often be eliminated through 
additional constraints on ^E". To enforce complete elim- 
ination at the design stage of the algorithm would re- 
quire non-local field constraints, which are difficult to 
implement directly in the optimization algorithm. For- 
tunately, for geometries involving only a few constraints 
(like the examples in Section III I, the optimized mag- 
netic potentials rarely contain any field zeros, as spuri- 
ous minima would reduce the achievable Zeeman pseudo- 
potential strength at the nominal trap positions. In more 
complex cases such as magnetic superlattices with mul- 
tiple traps per unit cell, or lattices of traps involving 
higher curvature constraints [for example from Eq. Q], 
it may become favorable for magnetic field zeros to oc- 
cur at points of special symmetry or even anywhere in 
the trap lattice. In these situations extra care should 
be taken in choosing appropriate additional constraints 
for the algorithm to eliminate or shift the spurious zeros 
from the trapping regions. 

It may be of interest to directly constrain the symme- 
try of the Zeeman pseudo-potential. As an example, the 
triangular lattice we present in Section pll A 2| has a high 
degree of symmetry with equal barriers between sites cho- 
sen to allow for symmetric tunneling. However, direct 
symmetry constraints on the potential are quadratic in 
the coefficients ma and thus cannot be incorporated into 
the linear form of the algorithm. This could be partly 
overcome through the use of nonlinear optimization algo- 
rithms which allow both linear and nonlinear constraints 
on the magnetic field and pseudo-potential respectively. 
We have applied a nonlinear optimization algorithm to 
some problems; however, we find it considerably slower 
and limited in the achievable resolution when compared 
with the linear programming algorithm. 



III. OPTIMIZED LATTICES 

In this section we present results for two-dimensional 
periodic lattices with a single microtrap per unit cell. As 



discussed in Section II B the simplest IP atom traps are 
constructed by constraining the magnetic field gradient 
to be of the form 



V (X R, 




1 I • 



(13) 



at the desired trap locations, where R^i^e^ip is a rotation 
matrix with three Euler angles specifying the orientation 
of the IP trap axis, ly — {sin 9 sin i/j, sin 6 cos ip, cos 9}, and 
the orientation (j) of the two perpendicular axes around 
this direction. The local field vector u can be left un- 
constrained since we are free to apply any external ho- 
mogeneous bias field Bq in order to null the local mag- 
netic field at the traps. We have run the optimization 
procedure for a large set of Euler angles and computed 
the corresponding optimal lattice geometries. In prac- 
tice we find that IP traps with in-plane loffe axis orien- 
tation (9 = 7r/2) provide good confinement, in contrast 
to out-of-plane traps {9 — 0) which typically result in 
weak or no confinement perpendicular to the surface [in 
the absence of curvature constraints, Eq. We stress 
however that the algorithm described in Section II C is 



general, and constraints on all derivatives of the scalar 



potential {ui, Vij 



and higher) can be applied to 



obtain any desired trap shape or orientation. In the fol- 
lowing we give a general and two specific examples of 
in-plane {9 — tt/2) IP magnetic lattices, chosen for their 
desirable properties for future planned experiments with 
ultracold atoms. 
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A. Two- wave lattices 

An interesting class of IP trap arrays with in-plane 
loffe axis can be described by the three-dimensional 
scalar magnetic potential 



cos a sm 



2ttx 



_27ry_ 
(itan(^ 

ci sin C 



(14) 



with a — tan 



^1 cos(-i/'+C) _ j|- fQ].j^g a^j]^ oblique lattice of IP 

cos ^ ^ 

traps above the points r — ni{c?, 0} -I- n2{(icosC,dsinC} 
with (711,712) G Z^, where d is the lattice period. Their 



loffe axis is {sin ip , cos ip , 0} , in accord with Eq. (13) for 
6 = tt/2 and cj) — 7r/4. The trapping height h is selected 
through the bias field, which must cancel the field 

, , , 27rh 2tt cos a r , 1 / s 

B(0,0,h) = -e"?-^ X — {-cosV',sinV',0} (15 

a COS?/' 

perpendicular to the loffe direction. What is special 
about these trap arrays is that there is a particular loffe 
field strength 



-e dsinC X 



TT COS a sin(2-0 + () 
cos^ Vcos(V' + C) 



(16) 



such that the Zeeman pseudo-potential ( 12 ) is invari- 



ant under exchange of the lattice axes (wallpaper group 
cmm), irrespective of ip. In particular this gives rise to 
equal tunneling barriers in the two lattice directions. It 
is interesting to note that this cmm symmetry is present 
neither in the magnetization pattern m{x, y) cx ^'(a;, y, 0) 
nor in the bias or loffe fields. Therefore, using the opti- 
mization algorithm of Section [II C| to find binary magne- 
tization patterns which optimally produce scalar poten- 
tials like Eq. (14 1, we can obtain experimentally relevant 



and non-trivial systems which are otherwise difficult to 
discover using manual methods. In what follows we per- 
form such calculations to construct approximately cmm- 
symmetric square and triangular lattices of IP traps. 



1. square lattice 

Square lattices are commonly produced using optical 
potentials but are a considerable challenge to implement 
with magnetic lattices, as the high degree of symmetry 
typically results in points of zero magnetic field strength 
which lead to Majorana spin flips [25 . However, from 
Eq. (14 1 with C, = tt/2 such square lattices of IP traps 
are readily generated with our algorithm. All values of 



ip except integer multiples of 7r/4 are allowed; angles 
close to 7r/4-|-ri7r/2, n € Z, give maximal trapping depth. 

A possible implementation of such a square lattice of 
IP traps is found by constraining the second derivative 




FIG. 1: Optimized magnetization pattern for a square lat- 
tice of IP traps (above blue points) constrained by Eq. (17 1. 
Trapping height is /i = d/2, and the loffe direction is given 
by the arrow [ip — 7r/3). The optimization yields C — 1.06 in 
Eq. (17 1. Blue (white) areas are fully magnetized (unmagne- 
tized) with m = 1 (m — 0). 



tensor of the scalar magnetic potential at the trapping 
sites to its value dictated by Eq. (14), 



C X 




cost/; 
— simp 




(17) 



The first and third derivatives u and w are left uncon- 
strained. In Fig. [1] we show 2x2 unit cells of the optimal 
magnetization pattern which generates this gradient ten- 
sor for Tp = 7r/3 and nominal trapping height h = d/2, 
i.e., half the inter-trap distance. The obtained pattern is 
typical in that it consist of large connected regions of uni- 
form magnetization with equal areas of magnetized and 
non-magnetic regions. Furthermore, the pattern resem- 
bles the 'staircase'-like structure recently experimentally 
demonstrated to create a two-dimensional lattice of IP 
traps [ini [3T]. This shows that the optimization algo- 
rithm can produce magnetization patterns which resem- 
ble manually optimized lattice configurations. 

The trapping height h' is chosen by varying the bias 
field strength and direction. For h! ^ h the curvature 



tensor may differ from Eq. ( 17 ); however, det t;(0, 0, h') = 
for any desired trap height, and therefore we are guar- 
anteed to be able to produce IP traps everywhere on a 
vertical loading trajectory. Details of a method for load- 
ing such trap arrays are described in Section [IV D[ 

Fig.[2]shows the associated Zeeman pseudo-potential in 
the trapping plane, brought close to cmm symmetric by 
using an loffe field strength which creates equal potential 



barriers in the two lattice directions [similar to Eq. ( 16 )] 



Notice that the cmm symmetry is not exact because t 
higher-order Fourier modes used to construct the binary 
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FIG. 2: Zeeman pseudo-potential [Eq.(12l] in the trapping 
plane {z = h = d/2) produced by the magnetization pattern 
of Fig. [l] at a lofle field strength which creates equal potential 
barriers (dashed separatrix) betwee n th e potential wells (blue 
points) in the two lattice directions ( 16 1. Darker colors signify 
deeper potentials. 



magnetization pattern (Fig. [ij perturb the simple ex- 
pression of Eq. (14); neither the magnetization pattern 
nor the loffe direction are cmm symmetric, and there- 



fore there is no a priori reason why perfect cmm symme- 
try should be achievable in the Zeeman pseudo-potential. 
However, for sufficient nominal trapping heights h/d this 
asymmetry may be negligible in practical applications (it 
is less than 3% of the potential range in Fig. [2]). Further, 
by a slight tuning of the loffe field strength we can force 
the effective tunneling rates for trapped atoms in the 
two lattice directions to be exactly matched. This lat- 
tice is highly tunable: by varying the external bias field 
and the loffe field strength the barriers in the x and y 
lattice directions can be tuned independently, allowing 
time-dependent control over anisotropic tunneling rates. 



2. triangular lattice 

Ultracold atoms in triangular lattices are of great ex- 
perimental interest as they exhibit a wide variety of novel 
quantum phases However the considerable chal- 

lenge involved in optically producing triangular lattices 
has so far limited experiments to a single recent realiza- 
tion [33]. The optimization algorithm is equally easily 
applied to produce triangular magnetic lattices, found 
from Eq. (14) with C = 7r/3. For triangular lattices it 
would be desirable to have a Zeeman pseudo-potential 




FIG. 3: Optimized magnetization pattern for a triangular lat- 
tice of IP traps (above blue points) constrained by Eq. 1 17 1 
with i/) = 57r/12 and additional field constraints at the poten- 
tial barrier given in the text. The nominal trapping height 
is /i = d/2, and the loffe direction is indicated by the arrow. 
The optimization yields C — 0.729. Blue (white) areas are 
fully magnetized (unmagnetized) with m = 1 (m = 0). 



but fully p^)m-symmetric in order to have three equiv- 
alent tunneling directions and a truly isotropic lattice. 
Unfortunately, as mentioned in Section [H C 1| such sym- 
metry requirements involve nonlinear constraints beyond 
the linear programming method. However, we find that 
sufficiently symmetric lattices can be obtained by apply- 
ing additional linear constraints at the barrier positions. 

In Fig. [3] we present a magnetization pattern which 
produces a useful triangular lattice of microtraps with 
equal barrier heights in three directions, and nominal 
trapping height h = d/2. The required constraints 
include the gradient tensor at the trap sites given by 
Eq. ( [l7| with ^ — 5n/12, as well as an additional field 
constraint at one of the barrier positions at {x, y, z} — 
{d/2, 0, d/2}. The required field constraint is 



',fioSM./d 



= C X (-0.0977), 



(18) 



which is not just cmm-symmetric, as given by Eq. (16), 



obtained automatically by nesting the linear program- 
ming algorithm within a nonlinear search algorithm aim- 
ing to null the pseudopotential barrier height difference 
(Fig. [4]). This two- level approach is possible due to the 
speed of the linear programming algorithm. Although 
the resulting pseudopotential is not fully p6m symmet- 
ric, the effective tunnelling rates between traps could 
be matched experimentally by varying the loffe field 
strength Bi and the trap height to independently tune 
the potential barriers in all three directions. 

As for the square lattice of the preceding section, the 
trapping height h' is chosen by varying the bias field 
strength and direction, and a vertical loading trajectory 
of IP traps can be produced (see Section IV D ) . 
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x/d 



Triangular lattice (Fig. [s]): A bias field of Bq = 
{2.9, 26.1, 0} G creates traps at /i = 2.5 /im, with 
Bi = 9.8 G. The barrier heights are equal in all 
three lattice directions, 34.6 G (2.32 mK) above the 
trap minimum. The trap depth away from the 
surface is 16.5 G (1.11 mK). The trap frequencies 
are = 27r x 150 kHz, ujy — 27r x 146 kHz, and 
uj^ = 2n X 36.9 kHz. 

These lattices provide sufficiently tight confinement to lo- 
calize the atoms to smaller than the optical wavelength. 
The Lamb-Dicke parameter is rjj = ^Wrecoii/i^j where 
Wfccoii = 27r X 3.771kHz for the ^'^Rb D2 line. For the 
transverse dimensions 77 < 0.18 and the atoms would en- 
ter the Lamb-Dicke regime {rf'{n) <^ 1) for mean vi- 
brational quanta {n) < 30. The magnetic fields and 
Zeeman pseudo-potential strengths for the given exam- 
ples scale with SAI^/d and the trap frequencies scale with 



FIG. 4: Zeeman pseudo-potential [Eq.(12l] in the trapping 
plane (z = h = d/2) produced by the magnetization pattern 
of Fig. |3] at a loffe field strength which creates equal poten- 
tial barriers (dashed separatrix) between the potential wells 
(blue points) in the three lattice directions ( 16 1. Darker colors 
signify deeper potentials. 



IV. EXPERIMENTAL CONSIDERATIONS 



A. Trapping parameters 



The lattice geometries of Section III provide high trap 
depths and strong confinement for magnetically trapped 
atoms. In the following we calculate the resulting trap 
parameters, assuming realistic parameters for a FePt thin 
film [T^ with Mz = 670kA/m and S = 0.3 /im (magne- 
tization current SM^ = 0.2 A) trapping ^^Rb atoms in 
the F — 2,mp = 2 hyperfine state of the 5^Si/2 elec- 
tronic ground state {gi = —0.501826 71(5) and 32 — 
0.499 836 43(5) in this level |35|). A small but optically 
resolvable lattice period of c? = 5 /im is chosen so that 
nearest-neighbor interactions such as the Rydberg block- 
ade effect would be observable [351 |3Zj ■ 

Square lattice (Fig. [T]): A bias field of Bq — 
{—45.6, 21.7, 0} G creates traps at ft, = 2.5 /fm 
above the film layer, with an loffe field of Bi — 
— 29.7G. The barrier heights are equal in both x 
and y directions, 37.9 G above the trap minimum. 
For *^Rb in the |2,2) state this corresponds to a 
barrier height of 2.55 mK. The trap depth away 
from the surface is 20.8 G (1.40 mK). Each trap 
is approximately cylindrically symmetric with long 
axis in the {1,1,0} direction. The trap frequencies 
are = 27r x 124 kHz, uj j_ = 2t: x 121kHz, and 
wii = 27r X 37.6 kHz. 



B. Microfabrication 

The proposed magnetic lattices can be easily produced 
using conventional microfabrication techniques. Success- 
ful methods used in the past have included magneto- 
optical recording [38U40j . hard-disk write head [20j . 
grooved substrates and uniform film coating 4L42j, laser 
ablation [33] , and optical or e-beam lithography followed 
by reactive ion etching [13 IS]- We favor lithography 
and ion etching since they provide greater freedom in the 
films used, arbitrary magnetization patterns are possible, 
and they can be readily extended to sub-micrometer res- 
olution. Lithographic resolutions < 50 nm are possible 
with e-beam lithography, but long exposure times may 
be required to produce reasonably large lattices. Opti- 
cal lithography is comparatively simple and can produce 
vast lattices; however the achievable resolution is limited 
by diffraction within the photo-resist layer to typically 
1 /im. 

The optimized magnetic patterns are highly compat- 
ible with lithographic patterning, as they involve large 
connected magnetized regions with smooth boundaries. 
The effect of finite resolution is qualitatively equivalent 
to truncating the Fourier series expansion of the magnetic 
scalar potential. Because the contributions of higher 
Fourier modes (corresponding to small features) on the 
magnetic pseudo-potential decay rapidly with distance 
from the surface [Eq. ([TT|)], the shape of the potential 
at the trap position is relatively insensitive to the fine 
details of the magnetization pattern. For a more quan- 
titative analysis, a perturbative study of the dependence 
of the pseudopotential on the magnetization in a given 
domain is readily performed with the parametrizations 
of Eqs. ([7]), (|9| etc. used in the optimization algorithm. 

We have compared the trapping potentials produced 
by the optimized lattices at full numerical resolution to 
calculations with truncated Fourier series. Taking the 
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triangular lattice as an example, we find that the lowest 
13 Fourier modes < 2 x ^) are sufficient to ad- 

equately reproduce the desired potential (Fig. This 
corresponds to an effective resolution of 7r/||fc|| = 0.25d. 
For the d — 5 /im period lattice this would require a res- 
olution of better than ~ 1.2 /xm, which can be achieved 
by optical lithography. With the reduced set of Fourier 
modes we find the pseudo-potential in the trapping plane 
maintains its symmetry and deviates from the full reso- 
lution calculation (Sec. IIIA2I by at most 1.5%. The 



trap depth increases by 0.4% while the barrier heights 
decrease by 0.2-0.4%; the trap frequencies increase by 
0.4%, 0.1%, 1.1% in the z, y, x directions. 



C. Effect of inhomogeneity 

The quality of the magnetic lattices may also be af- 
fected by film roughness or magnetization inhomogene- 
ity [171 m] ■ These lead to an additional randomly ori- 
ented magnetic field above the surface which causes small 
site-to-site variations in the potential experienced by the 
atoms. A simple model which takes into account random 
(white noise) fluctuations of the magnetization at a posi- 
tion above the edge of a magnetic film layer [JS] is used to 
estimate the potential energy variations between lattice 
sites. A typical rms surface roughness of around 5 nm 
and grain size of 50 nm would cause unwanted variations 
of less than 20 mG corresponding to energy variations of 
less than 30 kHz. A more significant contribution may be 
due to imperfect magnetization and canting of magnetic 
domains. For a typical remanent-to-saturation magne- 
tization ratio of 0.95 (squareness of the hysteresis loop) 
and a characteristic domain size of 50 nm, we anticipate 
additional variations of 60 mG (85 kHz). The estimated 
combined effect of surface roughness and magnetic inho- 
mogeneity is comparable to or less than the vibrational 
level spacing and is therefore not expected to be a se- 
rious obstacle for experiments. The anticipated site-to- 
site variations could be reduced further with high qual- 
ity Fe/Pt or Co/Pt multilayer thin-films, which can be 
atomically flat with small magnetic domains and higher 
squareness ratios [46] , 



D. Loading trajectory 

Atoms can be transferred to the proposed lattices sim- 
ply using an auxiliary macroscopic Z-shaped wire be- 
neath the chip surface [19 1 . An electric current through 
the Z-wire combined with an external bias field produces 
a single loffe-Pritchard magnetic trap [17], hundreds of 
micrometers above the film surface where the influence 
of the magnetic lattice is negligible. 

We have computed a loading trajectory which 
smoothly transfers atoms from the Z-wire trap to the 
magnetic lattice potential generated by the pattern of 
Fig. [3| and a lattice spacing of d = 5 fj,m. We assume 



a Z-wire positioned 0.4 mm beneath the magnetic film, 
with a central segment of L = 1 mm oriented such that 
its loffe axis matches that of the microtrap lattice at the 
nominal trap height h — 2.5 /im. In this symmetric situ- 
ation dot t>(0, 0, z) = Vz > for the total scalar mag- 
netic potential (lattice -I- Z-wire), and thus for any de- 
sired trapping height h' we can still find a homogeneous 
in-plane bias field Bq which creates IP traps t25j. By 
smoothly ramping down the Z-wire current or increasing 
the bias field strength, the trap minimum moves toward 
the surface and the Z-wire trap eventually merges with 
the lattice microtraps. 

For each trapping height between h' = 100 /zm and 
h' ^ h ~ 2.5 /im, the bias field (i?o) and Z-wire current 
{Iz) are optimized to maintain a trap depth equal to or 
greater than the final depth in the lattice of > 16.5 G 
to prevent loss of atoms during loading [Fig. [5ja)]. The 
loffe-field strength is chosen to start at Bi = 2 G and 
increases during the trajectory to its final value of Bi = 
9.8 G to create a symmetric triangular lattice (Fig. |4|. 
Initially the atoms are confined in the large potential well 
of the Z-wire trap at h! — 100 /im, however the super- 
imposed field of the lattice produces shallow secondary 
potential wells at z « 3/im [Fig. [5jb)]. The magnetic 
field strength as a function of height z during the loading 
trajectory is calculated at the position of the global mini- 
mum in the xy plane. As the fields are varied to move the 
Z-wire trap closer to the surface the minima approach, 
and at a distance of around /i' « 10 //m atoms can spill 
over to the secondary traps [Fig. [5](c)]. At this stage a 
thermal atom cloud would begin to fragment. Evapora- 
tive cooling in the combined Z-wire-|-lattice trap may be 
very efficient as the local phase-space density could be 
high in the small potential wells [48l - f50] . Finally the Z- 
wirc current is ramped off completely and the atoms are 
confined by the lattice potential alone [Fig. [sjjd)]. The 
precise number of traps that can be loaded in this way is 
difficult to estimate since it depends on the temperature 
of the clouds during loading and the rate of transfer. In 
practice, it seems possible to load a 1.5 x 0.5 mm" region 
of the lattice [21], corresponding to ~ 3 x lO"' populated 
microtraps. Throughout the transfer we ensure no field 
zeros are produced in the vicinity of the atoms, such that 
Majorana spin flip loss can be neglected. Additionally, 
the field strength at the film surface is ||i?|| > 120 G in 
the z = \ fxni plane, corresponding to a repulsive atom- 
surface barrier much greater than the attractive Casimir- 
Polder potential, thus preventing atom loss to the sur- 
face. 



V. CONCLUSIONS 

We have introduced a linear programming algorithm 
tailored to the problem of designing two-dimensional 
magnetic lattices of loffe-Pritchard traps for ultracold 
atoms. Previously, the design of such lattices required 
substantial experience and trial-and-error. The algo- 
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difficult to obtain using manual methods. 

The strengths of the algorithm were exemplified in a 
desirable square magnetic lattice and in the new case of 
a triangular magnetic lattice potential. The generated 
magnetic patterns consist of smooth connected regions 
of uniform magnetization, which are easily fabricated us- 
ing existing methods. With realistic magnetic-film pa- 
rameters these lattices will produce vast arrays of mi- 
croscopic IP traps with tight confinement and high trap 
depths. The resulting potentials are of great interest for 
studies of quantum gases in lattice potentials [3 , 
condensed-matter analog systems [H| , or quantum 

information processing with neutral atoms [9Vll4j . 

The similarities between the algorithms for design- 
ing magnetic loffe-Pritchard traps and electric radio- 
frequency (rf) ion traps could be exploited to produce 
hybrid traps, as a natural route to confine individually 
controlled rf-trapped ions within magnetically trapped 
ultracold quantum gases [5T1[S2]. By appropriate changes 
of the external bias field the neutral atoms could then be 
moved to overlap with the ions or separated from them 
in a well-controlled manner. Further, once the ions are 
cooled to their motional ground state the rf traps can be 
switched off, transferring the ions to the magnetic mi- 
crotraps (provided that their magnetic moments are ap- 
propriate) in order to eliminate rf micromotion [531 154) . 
In practice one possible approach would be to add extra 
constraints in Eq. ( 10 1 to produce patterns that simul- 
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taneously provide magnetic traps for neutral atoms and, 
by applying rf voltages directly to the proposed magnetic 
structures (provided their rf impedance is suitable), rf 
traps for ions at the same sites or at sites of our choice. 
Alternatively one can pattern an appropriately shaped 
conductive rf layer on top of the magnetic structures, 
thus decoupling the patterns for atom and ion traps. 



FIG. 5: Calculated trajectory for loading atoms to the tri- 
angular lattice, (a) Z-wire current 1^ (solid line) and ap- 
plied bias field Bo (dashed, dotted, and dash-dotted lines 
for the X, y, and z components, respectively) as functions 
of trapping height h' . Magnetic field strength ||B|| as a func- 
tion of distance from the film surface: (b) prior to transfer 
{h' — 100 /xm), (c) during transfer {h' = 10 /xm) and (d) after 
transfer to the magnetic lattice is complete {h' = h = 2.5 ^m). 
Note the logarithmic z scale. 

rithm automatically generates single-layer binary mag- 
netization patterns which produce desired lattice sym- 
metries with specified trap parameters. It allows for de- 
signing non-trivial geometries which would be extremely 
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